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I. INTRODUCTION 



Recently the ppssibilitv to measure ^13 by a reactor experiment has attracted much at- 
tention B 3 3 B II B 0, i, 0. To achieve sensitivity sin^ 2^13 ~ 0.01, reduction of the 
systematic errors is crucial, and near and far detectors seem to be necessary for that purpose. 
On the other hand, it appears to be advantageous to do an experiment at a multi-reactor 
site to gain statistics and high signal to noise ratio, and in fact in the most of cases consid- 
ered in jH 3 3 3 3 3 Hi) there are more than one reactor. In this paper we discuss 
the systematic errors in reactor neutrino oscillation experiments with multi reactors and 
multi detectors in an analytical way. In Sect. 2 we discuss the cases with a single reactor 
to illustrate our analytical method. In Sect. 3 we consider the cases with reactors and 
show that the larger Ur gives totally the smaller contribution to the sensitivity from the 
uncorrelated errors of the fluxes. Irrespective of the number of reactors, if there are more 
than one detectors, we can cancel the correlated errors which includes the error of the fluxes. 
We emphasize in this paper that the sensitivity on sin^ 2^13 with vanishing statistical errors 
(we refer to the sensitivity as the systematic limit) is dominated by the uncorrelated error 
of the detectors in most cases. It is also emphasized that a lot of caution has to be exer- 
cised to estimate the uncorrelated error. In the appendix we give some details on how to 
derive the analytic results used in the main text, using the equivalence of the pull method 
and the covariance matrix approach jiol . 11,0, 3]- Throughout this paper we do not use 



the binning of the numbers of events because the discussions on the uncorrelated bin-to-bin 
systematic errors are complicated. Also we will discuss only the systematic errors, i.e., we 
will consider the case where the statistical errors are negligibly small. 



II. SYSTEMATIC ERRORS 

To discuss the systematic limit on sin^ 26'i3 iii neutrino oscillation experiments with mul- 
tiple reactors, we have to introduce the systematic errors of the detectors and the reactors 
(fluxes). There are two kinds of systematic errors of the detectors, namely, the correlated 
error a^. and uncorrelated error cr„. The former includes the theoretical uncertainty in the 
cross section of detection, etc. while the latter does the uncertainty in the baseline lengths, 
a portion of measuring the detector volume, a part of the detection efficiency, etc. As for 
the systematic errors of the reactors (fluxes), the correlated error Uc^^ consists of the uncer- 
tainties in the spectrum of the flux, etc. whereas the uncorrelated error au consists of 
the uncertainties in the composition of the fuel, etc. In this paper we adopt the reference 
values for cT p an d used in B, where basically the same reference values as in the Bugey 
experiment |l^ were assumed. cXc and can be estimated to be 

cTu = 0.8%/V2 = 0.6% 

a, = ^ (2.7%)2 - (2.1%)2 - (0.8%/V2)2 = 1.6%, (1) 

where the factor appears because the relative normalization crrei=0.8% in [l3] is related 
to (Tu by (Jrei = ^Oy,. In the estimation of cTc, we used 2.7% total error and 2.1% error of the 
flux which are the values in the CHOOZ experiment. As for the correlated and uncorrelated 
errors of the the flux from the reactors, we adopt the same reference values as those used 



2 



by the KamLAND experiment [l5| 

a« = 2.5%, 

a« = 2.3%. (2) 

Note that the word "correlated" means just the type of the error, and then correlated errors 
exist even if there is no partner. 



III. ONE REACTOR 

To explain our analytical approach, let us start with the simplest case, namely the case 
with one reactor. 



A. One detector 



Let m be the measured number of events at the detector, t be the theoretical predictions 
(hypothesis) to be tested, is defined as 



X = mm 

a's 



m 



t(l + ac + ac^ + On 



-\ 2 



tCTn 




0-, 



(3) 



where ttc, c^c^ and au^ are the variables of noises to introduce the systematic errors (T„, 



(Jc, ac and au , respectively. We give an easier way to derive (jHI) in the Appendix 13 where 
integration over the a variables as those of Gaussian, instead of minimizing with respect to 
these variables, do the same job. Eq. Q shows that the square of the total systematic error 
is given by the sum of the squares of all the systematic errors. 

Our strategy in this paper is to assume no neutrino oscillation for the theoretical pre- 
dictions t's and assume the number of events with oscillations for the measured values m's. 
Then, we examine whether a hypothesis with no oscillation is excluded or not, say at the 
90%CL, from the value of x^- In the context of neutrino oscillation experiments, we have 



m 
T 



sin 26 ( sin 



4:E 



(4) 



in the two flavor framework, where 6 is the mixing angle, Am^ is the mass squared difference^, 
E is the neutrino energy, L is the distance of the reactor and the detector, and 

2 / Am^L/ 

V 4:E 



sm 



dE e{E)f{E)a{E) sin^ 



dE e{E)f{E)a{E) 



^ Throughout this paper, we use the two flavor framework. To translate it into the three flavor notation, 
and Am^ should be interpreted as ^13 and |Am|]^|, respectively. 
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e{E), f{E), a{E) stand for the detection efficiency, the neutrino flux, and the cross section, 
respectively. 



B. Two detectors 



Next let us discuss a less trivial example with a single reactor, one near and one far 
detectors. Let and mf be the measured numbers of events at the near and far detectors, 
tn and tf be the theoretical predictions, respectively. Then is given by 



X 



u 




1 2 



+ 



iTLf — tf(l + etc + ai'^ + an 



(5) 



where we have assumed that the uncorrelated errors for the two detectors are the same and 
are equal to a^- Eq. ^ can be evaluated also by integrating over the variables Oc, etc. as 
Gaussian instead of minimizing with respect to these variables. After some calculations (See 
Appendix 1X1 for details), we obtain 



X 



tn 



1 



-1 



where 



V 



W^2 



0",, 



or, 



Wa2 



(6) 



is the covariance matrix; I2 represents 2x2 identity matrix and H2 does a 2 x 2 matrix 
whose elements are all unity. It is seen that only the covariant matrix in the depends 
on the errors. Note that any liner transformation of V does not change the value of x^- 
Diagonalization of V is, however, worthwhile to investigate analytically the behavior of x^- 
After diagonalizing V we have 



X 



[{mJU - 1) + Wtf - 1)]' , [{mJU - 1) - {rrif/U - 1)] 



4a2 + 4(cri''^)2 + 4((T, 



2al 



+ 



2al 



(7) 



sin^ 29 



{D{L,) + D{L^)y 



+ 



[D{L,) - D{K)f 
2ol 



(8) 



_4(T2 + 4K0' + 4K0' + 2a2 

where and Lf are the distances from the reactor to the near and far detector, respectively; 
We have defined 



D{V) 



\^sm 



4E 



(9) 
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The first term on the right hand side in Eq. ((Tj) stands for the contribution from the sum of 
the yields at the near and far detectors, while the second term corresponds to the difference 
between them. The first term determines the normalization of flux, namely the sensitivity 
on sin^ 2^ at very large |Am^| where all D{L) becomes 0.5sin^2^. On the other hand, the 
second term gives the main sensitivity at the concerned value of |Am^| (e.g. 2.5 x 10~'^eV^) 
as we see below. 

Putting the reference values ((H) and Q together, we have 



2ai = (0.8%) 
4ae^ + 4(aW)2 + 4(a«)2 + 2a^ = (7.6%) 
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We can ignore the contribution from (4(T^ + A{a-u^Y + A{ac^)'^ + 2ct^) ^ in Eq. (jHJ because 
that is only 1% compared to that from {2a'^)^^.^ Hence, is given approximately by 

,^.,i„4 29(£(M^£(M!. (10) 

We see that the main sensitivity is determined indeed by the relative normalization error 
(Jrei = v^o"u. By comparing Q and (HU)), it is clear that the sensitivity is improved sig- 
nificantly by virtue of near detector. The hypothesis of no oscillation is excluded at the 
90%CL if is larger than 2.7, which corresponds to the value at the 90%CL for one degree 
of freedom. This implies that the systematic limit on sin^ 29 at the 90%CL, namely the 
sensitivity in the limit of infinite statistics, is given by 

(sin^ 2^^''' a/27 (11) 

1^1^ ^'^^ limit n(TA _ n(T \- ^^^^ 



D{U)-D{L^y 

Eq. (fTTj) also tells us that, in order to optimize (sin^ '^^Tvixait^^ ^ given value of (Tu, we have 
to maximize D{L{) = (sin^ (Am^Lf/4£')) while minimizing D{L^) = (sin^ (Am^Ln/4i?)). 
Note that D{Li) can not be unity because of neutrino energy spectrum; The possible maxi- 
mum value of D{L{)—D{Ln) is 0.82, which is attained for Am? = 2.5 x 10~^eV^, Lf = 1.8km, 

and Ln = 0. Then, we can estimate the lower bound of (sin^ 2^) at a single reactor 

experiment, assuming that the uncorrelated error is smaller enough than other errors: 

In practice, however, D{L^) will not be able to vanish. Assuming Am^ = 2.5 x 10^'^eV^, 
Lf = 1.7km, and = 0.3km, we have D{L{) ~ 0.82 and D{Ln) ~ 0.07. Then, = 
(0.8/\/2)% gives the sensitivity 

(sin2 2e)7^°°'^ ^ VT7-— ^^^^^ ^ 3.1a, ^ 0.018. 
This sensitivity corresponds to and agrees with the value obtained numerically in 



^ This is more or less the derivation of used in 



5 



IV. Ur REACTORS 



It is straightforward to generalize the argument in the previous section to a general case 
with multi-reactors and multi-detectors. The covariance matrix V is given by (T„x(unit 
matrix) + (the rest), and in most cases, as long as the near detectors are placed properly, the 
determinant of (the rest) is zero or very small compared to a^. The minimum eigenvalue of 
the covariance matrix, which gives main contribution to x^i is approximately given by a^. 
Therefore, the systematic limit is dominated by the uncorrelated error Uu also in general 
cases. 



A. One detector 



As in Sect. IIIIl as a warming up, let us consider the case with one detector and multiple 
reactors (See Fig^a)). When there are nr{> 1) reactors, the total number m of the mea- 
sured events is a sum of contributions nia (a = 1, ■ ■ ■ , rir) from each reactor, and this is also 
the case for the theoretical predictions t and (a = 1, ■ ■ ■ , rir). So we have 



m 



a=l 



a=l 



Taking these systematic errors into consideration, the total number t of the theoretical 
prediction is redefined as 



a=l 



t (1 + + «w) + ta(y^^ = t 1 + «c + a^:^ + 5Z f « 



(r) 



a=l 



a=l 



where aud is the variable to introduce the uncorrelated of the flux from the a-th reactor, 
is defined as 



X = min 

a's I r^O"^ 



m 



a=l 



2 / (r) \ 2 / (r) ^ ^ 



cr, 



a=l \Cru 



where we have again assumed for simplicity that the size of the uncorrelated error in the 

(r) (t) 

flux from the reactors is common: ana = (Tu ■ Performing minimization with respect to the 
a variables, we get (See Appendix 



X 



m 

1 

t 



a=l 



(13) 
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By comparing ()13|) with Q, we find that the multiple reactor nature affects (the sensi- 
tivity on sin^ 2^) through a^^. Since there is only one detector in this case, the correlated 
error was not canceled in aiid it contributes to the systematic limit on sin^ 26. However, 
if the yield from each reactor is approximately equal, i.e., if 

7 - (14) 

t Ur 

then the contribution of the uncorrelated error of the reactors becomes 



a. 



Comparing Eqs. and (jl3j) . we observe that the contribution of the uncorrelated error of 
the reactors decreases as the number of the reactors increases, as long as the condition (jl5p 
is satisfied.^ This is because the average of independent fluctuations is smaller than a 
single fluctuation. 



B. detectors 



Let us now discuss more general cases with reactors and detectors. For simplicity 
we assume again that the size of the uncorrelated errors for the detectors are the same and 
the size of the uncorrelated errors in the flux from the reactors are also the same: a^j = o"u, 
(Tua = cr". Let tja {mja) bc the theoretical prediction (measured value) for the number of 
events of neutrinos from the a-th reactor (a = 1, ■ ■ ■ , n^) at the j-th detector (j = 1, ■ • ■ , n^) 
and tj = J2a=i ^ja ("^i = ^^=1 "^ja) be the theoretical (measured) total number of events 
at the j-th detector. Then generalizing the discussions in the previous sections, we have 

2 



X 



mm 

a's 



rid 



1 ^i*^" 



rrij - tj 



a=l 



a 



(r) 



+ I — 



+ 



a 



a=l 




After some calculation (See Appendix [X]) , we have 



X 



mi _ 



m. 



rid 



''nd 



( ^-1 \ 



m. 



nd 



nd 



where 



a=l 



tj tfc 



(16) 



^ One can show from the Cauchy-Schwarz inequahty that J22=i i^a/t) < 1 always holds even if the con- 
dition H15|l is not satisfied. Hence, the contribution of the uncorrelated error of the reactors decreases 
always. 
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1. rir reactors and {n^ + 1) detectors 



In Sect, mil we have seen that the correlated error is canceled in the case of a single 
reactor experiment with one near and one far detectors. Now we would like to ask the 
following question: what happens to this cancellation in the case of an experiment with 
multiple reactors and detectors? To answer this question, let us consider the ideal case with 
rir reactors and (n^ + l) detectors, where each reactor has a near detector in its neighborhood 
and each reactor produces the same number of events at a far detector (See Figmb)): 



near detectors : 
far detector : 



t 



ailr + l 



1 

rir 



In this case the element of the covariance matrix becomes 

/ 1 

■•• 



V = afj, 



U ^Ur + l 



+ 







lln.,\ 





1 



I /rir ) 



(17) 



(18) 



where /«,.+! is an [ur + 1) x {rir + 1) unit matrix, and -ffn^+i is an (n^ + 1) x (jir + 1) matrix 
defined by 



rir+l 



(19) 



Here we assume the following conditions^: 



4E 



«1, 



AE 



TC 



(20) 



These conditions have to be satisfied in an experiment which aims to measure ^13. From 
Eqs. (1201) we have 



/ mi/ti - 1 \ 







- sin^ 29 ( sin^ 



AE 



sm'^2eD{Li 





Vi/ 



^ For simplicity we assume here that the distance between the a-th reactor and its near detector is equal to 
Ln for a = 1, • • • ,nr- In order for H17|l to be satisfied, Eq. (|2()|l is necessary. So in this ideal situation which 
we are considering, the dependence on (sin^ (Am^Ln/4£')) cannot be discussed in a manner consistent 
with the assumption l{T7|l. 
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and becomes (See appendix iBl for derivation.) 



X 



sin^ 2eD{U 



rir + 1 



al + {rir + 1) 



(r) 



^(0 



+ 



rir + l al 



.(21) 



Then, the systematic hmit on sin^ 26 at the 90%CL is given by 



(sin 2^)i.^.^ 




1 

'2 



(r) 



.(0 




(22) 



(23) 



For example, we obtain (sin^ 29y^^^^^ ~ 0.012 by assuming seven reactors, L{ = 1.7km, and 

(Tu = 0.8/ \/2%. As in the case with one reactor, the dominant contribution to the systematic 
hmit comes from the uncorrelated error a^. The contribution of the uncorrelated error of 
flux, au \ is reduced in by a factor of due to the averaging over the independent rir 
fluctuations; Although this reduction is a potential merit of the multiple reactor complex, 
it is irrelevant to the sensitivity because such an effect comes in the correlated error among 
detectors which is almost canceled in the multi-detector system. The factor ^^/T^^TJn^ 
which appears in the dominant contribution by indicates that the effective systematic 
error decreases as the number (1 + rir) of the detectors increases, since more information is 
obtained with more detectors. To conclude, the answer to the question at the beginning of 
this subsection is that the cancellation of the correlated error occurs also in the ideal case 
with rir reactors and {ur + 1) detectors, and once again the systematic hmit is dominated by 
(Tu. It should be noted that the number rir of the near detectors in this case is sufficient but 
not necessary to guarantee this cancellation, as we will see below in the case of the KASKA 
plan. 



2. The case of the KASKA plan 

The Kashiwazaki-Kariwa nuclear power station consists of two clusters of reactors, and 
one cluster consists of four reactors while the other consists of three (See FigHl^c)). According 
to the discussion in the previous section, we understand that near-far cancellation occurs for 
the KASKA case if we have seven near detectors. In the KASKA plan, however, not each 
reactors but each cluster of reactors assumed to have a near detector. In this subsection 
we would like to clarify the following questions on the KASKA plan Q: (a) Is the number 
of near detectors sufficient for the cancellation of the correlated error? (b) What are the 
effects of multiple sources? (c) Is the KASKA plan optimized with respect to the sensitivity 
to sin^ 26'? Here we again assume that the size of the uncorrelated error in the flux from 
the reactors is common and the size of the uncorrelated errors of the three detectors are the 
same. 
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Before we discuss the systematic limit for the actual KASKA plan, let us consider the 
ideal limit, in which all the reactors in each cluster shrinks to one point as is shown in 
Fig^d); The ideal limit is similar to the case discussed in the section llVB ![ namely the 
case of two reactors and three detectors. In this ideal limit, we have 



near 1: 



near 2: 



ta2 



I,-- - ,4) 

= 5,6,7) 
I,--- ,4) 

= 5,6,7) 



(24) 



(25) 



far: 



ta3 



and the covariance matrix is given by 



V 



(r)^2 



.7) 



1/4 1/7' 
1/3 1/7 
1/7 1/7 1/7 



0"; 



v 



«)^ 



+ 



(r)N2 



al + (a«)^ + ^ 



Wa2 



3 

Wn2 



7 



+ (a«)^ + 

Wa2 



a. 



[(Jc 



+ 



al + (a«)2 + 



Wn2 



7 

(r)A2 



0- 



7 

Wn2 



(26) 



(27) 



where J3 is a 3 x 3 unit matrix and ifs is a 3 x 3 matrix defined in Eq. ()19|1 . After diagonalizing 
V ^ we obtain the systematic limit (See appendi^O) 



(sin 2^),^. 



sys only 
limit 



(28) 



" D{Li) 7 

o"u is the dominant contribution to the systematic limit in Eq. (j28p and the correlated errors 
are canceled due to the near-far detector complex. The reason that we have the factor 
1^74/49 instead of (1 + l/A^)^^^|Af=2 = a/3/2 is because the ratio of the Ve yield at the first 
cluster to that at the second one is 4:3 instead of 1:1 assumed in (fT7|) . 

In reality, however ^the conditions P^ - ()26|) are not exactly satisfied in the setting of the 
actual KASKA plan jj]. Let us evaluate the exact eigenvalues of the covariance matrix by 
taking into account the actual parameters in 3. Table HI shows the power of the reactors and 
the distance between the seven reactors and the three detectors. From this we can calculate 



the fraction taj/tj (a = 
matrix is now given by 

V = alh + 



7, j = 1,2,3) which is given in Table UTl The covariance 



0.221 0.042 0.149 
+ (^c'T ^3+ W"))' I 0.042 0.298 0.138 

0.149 0.138 0.145 



(29) 
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Diagonalization of V can be done only numerically with the reference values in Eqs. (^J and 
(j21). We find that the minimum eigenvalue of is 1.04 x cr^. The eigenvalue shows that 
the cancellation of the correlated errors occurs also in the actual KASKA plan though the 
number of near detectors is less than that of reactors. The systematic limit on sin^ 29 is 
approximately given by the contribution from the minimum eigenvalue (See appendi^Hl): 



(sin^ 2^)^.°"^^ Vl^^^ ^ 3.9^^ ^ 0.022, (30) 

V /limit „ . , ^ ni „ . , X ni „ . ^ X ^^ ^ ' ^ ^ 



where u^^^ = {ui\u'2\u^^'*)'^ is the eigenvector of V which corresponds to the minimum 



eigenvalue (~ cr^) of V , and D{Lj) is the average of each contribution D{Laj): 

^ = E ^ (^i^' (^IF^) ) = 3(fer)). (31) 

Here Laj is the distance between the a-th reactor and the j-th detector, and taj/tj is the 
fraction of the yield from the a-th reactor at the detector j=l, 2, 3. When the near-far 

cancellation occurs sufficiently, the value of (sin^ 26')j^.°°''^ /cTu gives a good measure for 
the power of a reactor experiment almost independently of assumptions of error sizes; The 
smaller value means the better setup of reactor experiments. 

To see how effectively the correlated errors are canceled in the actual KASKA plan, 
comparison is given in Fig|21 between the sensitivities to sin^ 29 of the actual KASKA plan 
(Fig^c)) and of a hypothetical experiment with a single reactor and two detectors (300m 
and 1.3km baselines) depicted in Fig^e)); The same value of uncorrelated systematic error 
au and the same data size (=20t-yr) are used for each case. We observe that there is little 
difference between the sensitivities at Am^ = 2.5 x 10~^eV^. Also it is remarkable that the 
sensitivity of the actual KASKA plan for higher value of Am^ is better than of the single 
reactor experiment. This is exactly because of the reduction of the uncorrelated error from 
due to the nature of multi reactors (cf. Eq. (fT3j) ). where the near detectors play a role as far 
detectors in this case. Here, it should be mentioned that we see in Fig. Elthat the sensitivity 
in KASKA changes only to sin^ 29 ~ 0.03 even for Am^ = 2 x 10"^eV^. 

Now let us see the effects of multiple sources on the systematic limit (sin^ 26')j'^.°°''^ in 
the KASKA plan. Since the effect comes with au^ (compare (fT^ with ©), we can ignore 
other errors. The systematic limit for the case is presented in Fig. |21 with solid line; cxu 
is assumed to be 2.3%. Note that the systematic limit is extremely close to zero for one 
reactor case because of ideal near-far cancellation. The solid line in Fig. |21 shows, however, 
that about 0.4% error remains. Thus, we find that the effect of the nature of multiple 
sources is only about 0.4%. It is vary small compared to the relative normalization error 
\/2(ju(~ 0.8%) of the detectors.^ Once we know that two near detectors are sufficient 
for near-far cancellation with the Kashiwazaki-Kariwa nuclear power station, we should 
investigate the optimal locations of the near detectors for the cancellation. While the two 
near detectors are assumed to be very close to the reactors in the ideal limit (Fig. ^d)), 
the near detectors in the actual KASKA case cannot be too close to the reactors. If each of 



^ Note that different types of errors are combined as the sum of squared: 0.8^ + 0.4^ ~ 0.89^. 
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the two near detectors is too close to a reactor, it is impossible to cancel the uncorrelated 
error of the flux from other reactors, namely an'' ■ Hence, the optimization of the locations 
of two near detectors at the Kashiwazaki-Kariwa site is nontrivial. At first, in order to see 
the importance of the location of near detectors, we compare the systematic limits of the 
actual KASKA plan (Fig. mc)) with that of a hypothetical experiment (Fig. HJ^f)) in Fig. El 
In the hypothetical experiment, one near detector is very close to the reactor 7^1 in the first 
cluster while the other detector is very close to the reactor ^5 in the second one. In Fig. El 
we use (Tu = 2.3% and all other systematic errors are set to zero because cxu ■* dominates the 
difference of near-far cancellations in the two cases. For Am? = 2.5 x 10~^eV^, we obtain 
(sin2 2^)'f^.°°'^ ~ 4 X 10^3 (2 x lO^^) for the actual KASKA plan (the hypothetical case). 

(t) 

In the hypothetical case, the sensitivity is deteriorated because the correlated error au is 
not canceled sufficiently and about 2% error remains. 

We can extend this analysis to that with the arbitrary position of each detector. To do 
the analysis, we first obtain the optimized positions of the detectors. And then we examine 
the sensitivity to sin^ 26 by varying the position of each detector, leaving the locations of 
the remaining detectors in the optimized ones. The results are given by the contour plots 
in Fig. m without statistical errors and in Fig. El with the data size of 20 ton-yr, where the 
reference values in Eqs. (|H) and Q are used. In these figures the locations of the detectors 
are also depicted for the optimized case and for the currently planned case. From these two 
figures we observe that the distance between each near detector and the reactors in each 
cluster is approximately (300±130)m in the optimized case. This results in slightly poorer 
sensitivity to sin^ 29, compared with the hypothetical single reactor case with a near detector 
which is arbitrarily close to the reactor (See Eq. (fT^ with = 0). From Figs. El IH El we 
see that the positions of the near detectors in the KASKA plan are appropriate for the near- 
far cancellation and almost optimized. Therefore, it does not suffer from the multi-reactor 
nature such as the variety of reactor powers or distances between the reactors and the near 
detectors. 

V. DISCUSSION AND CONCLUSION 

Using the analytical method, we estimated the systematic limits (sensitivity without sta- 
tistical error) on the neutrino oscillation parameter sin^ 29i3 in various setups of reactor 
experiments. In the simplest case, where there is one reactor and two detectors, the corre- 
lated systematic error is canceled. In the case of multiple rir reactors, we showed that the 
correlated systematic errors of the reactors and the detectors are canceled as naive expecta- 
tion if the number of detectors ria is sufficient [rid = rir + l). We found that multiple reactors 
and detectors set up has an advantage of the reduction of the remaining uncorrelated error 
if the set up is appropriate for the near-far cancellation of correlated errors. On the other 
hand, we explicitly showed that the contribution to the sensitivity to sin^ 26'i3 from the the 
uncorrelated error of the flux, which controls the multi-reactor nature, is negligibly small al- 
though there are only three detectors for seven reactors {ud < rir + l). The only disadvantage 
of experiments with < + 1 is that one cannot put the near detectors arbitrarily close 
to one of the reactors (even if one neglects the technical difficulties), because that would 
ruin the cancellation of the uncorrelated error of the flux, as we have seen explicitly in the 
KASKA case. We presented also the optimal positions of detectors in the KASKA plan; The 
planned position of near detectors are close to the optimal ones although it is better if the 
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baseline length for the far detector becomes longer beyond the bound of the power station 
site. In all cases studies here, it is the factor (j^j [{|sin^ {^m^Li/AE)^ — (svo? (Am^Ln/4ii^))] 
that sets the systematic limit on sin^2^i3, and hence it is quite important to estimate (Ju 
carefully. The factor (sin^ '^^T\im\t^^ 1^^ seems to be a good measure of the power of a reactor 
experiment almost independently of assumptions of error sizes. 



APPENDIX A: DERIVATION OF THE COVARIANCE MATRIX 

In this appendix we first show that the form of which is expressed as the minimum of 
the function of the a variables with respect to these variables leads to the form of which 
is bilinear in the ratio Trij/tj of the measured value rrtj_ divided by the theoretical prediction 
tj. This has been known in the literature [13, U, E^, as the equivalence between the 
so-called pull approach and the covariance matrix approach. And then we show that the 
same job can be done by integration of exp(— over the variables rrij/tj. 

In the cases which we are considering, the correlated systematic errors t\al^i, ■ ■ ■ ,t^cr^£ 
are introduced by the variables 



a 



Introducing the notation 



Vj 




can be written as 



X 



mm 



{y-Hay D-^ {y - Ha) + a' D;'a 

{a - A-^HD-^yf A{a - A-^HD^^y) + f {D'^ - D'^ H A'^ D'^) y 
f {Dn' - D-'HA-'H^D-^) (Al) 



mm 

d 



where 



H 



is an n X £ matrix. 



= diag(a^i,-- - 



is an n X n diagonal matrix whose element is the normalized uncorrelated systematic error 
CTyj- for the variable (j = 1, ■ ■ ■ , n), 

Dc = diag(a,\,--- ,al^) 
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is an i X i diagonal matrix whose element is the normalized correlated systematic error cr^^- 
for the variable aj {j = 1, ■ ■ ■ ,i), and we have defined 



A = D-^ + H^D-^H. 

Note that we can incorporate the effect of the statistical errors in our formalism by redefining 
Du — > diag((Ty]^ + 1/y/ti, ■ ■ ■ ,<jIi + l/^/t^), although we do not discuss the statistical errors 
in the present paper. From Eq. ()A1|) we see that the covariance matrix V is given by 

V={D-'-D-'HA-'H^D-y\ 

We could prove by brute force that V can be written as 

V = D^ + HD,H^, (A2) 

but it is much easier to prove it by expressing the matrix element Vij as the integral of 
exp(— over the variables yj. 

First of all, let us prove that the matrix element V^j can be written as 

Vij =-^y J dyViVj exp (^-^y'^V'^i/^ , (A3) 
where N'y is the normalization constant defined by 

Proof of Eq. ()A3j) goes as follows. Diagonalizing the covariance matrix V, which is real 
symmetric, by an orthogonal matrix O 

V = O^DO = O^diag ■ ■ ■ , Vn) O, 

the exponent can be rewritten as 

y^V-'y = y^O^D-'Oy = y''^D-'y', 



so that we have 



A/'y j dyyiyj exp (^-^V^V 
(O),, {D\, {0\^ 



Vr 



Thus Eq. ()A3|) is proved. 
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Now Eq. ()A3|) can be simplified by expressing as the integral over the variables a of the 
original x^'- 



MyUa I dy I ddyiV, exp <j [(a - A-^HD~^yf A{a - A-^HD-'y) 



+f {D-' - D-'HA-'H^D-') y] } 



da y^yj exp <! {y - Haf D^^ {y - Ha) + a^D^ ^a \ , (A4) 



where we have used Eq. ()A1|1 . the normalization constant TVq, is defined by 

1 



da exp < — 



2 L 



{a - A-^HD-^yf A [a - A-'HD-'y) 



and Afy and Afa are related by 



AfyK = (27r) (det D^f (det D,f^ . 

Eq. ()A4|) can be easily calculated by shifting the variable y ^ y" = y — Ha: 



Vij = AfyAfo, j dy'' j da (y' - Ha^ (y" - Haj 

exp < — - {y — Ha) D~ {y — Ha) + a D~ a 

= m., + (H),, (/^c)h 

= {D^ + HD,H^)^^. 

Hence Eq. ()A2|) is proved. 

In the case of one reactor with one detector (cf. Eq. ©), we have 

= diag {al 
H = (1,1,1), 

so the covariance matrix is 

V = D^ + HDM^ = al + al + (aW)^ + 
In the case of one reactor with two detectors (cf. Eq. ©), we have 



^ " " ' 1 



= diag {al (a«)^ (a«)^) 

1 1 1 
1 1 1 



H 
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so that we obtain 
V = + HD^H^ 



In the case of reactors with one detector (cf. Eq. (fT^ ). we have 



H 

so that we obtain 

V = D^ + HD,H^ 



= diag [al (a^) 



< + 1 1, 1, - 



ti tr, 



f 1 \ 

1 

ti/t 

\ tnjt J 



In the case of Ur reactors with Ud detector (cf. Eq. (HHI)), we have 



= diag [al 



/ 1 1 *L ilill \ 

' i i ' 



1 1 ^"r 
t ' ' ' t 



where . is an x unit matrix, so that we obtain 



/l 1 ^ 



t ' ' ' t 



\ 



a=l 



( 1 
1 



diag[a,^(a«)^(a«)^■■■,(a«f] 



where if^i^ is an x ra^ matrix whose elements are all 1 (cf. Eq. ()19j] ). 



1 \ 

1 



tilt ■■■ ti/t 

\tnjt ■■■ tnjtj 
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APPENDIX B: DERIVATION OF (EU) IN THE CASE WITH REACTORS 
AND Ur + 1 DETECTORS 



V in Eq. (fT^ can be diagonalized as 

U~^VU 



diag < {rir + 1) 



2 I (Tu 



where t/ is a unitary matrix defined by 



U = {u^'\ 



u 



(1) 



1 



(2) 



1 

7! 







u 



<3) 



1 

7! 



1 

-2 




u 



{rir+l) 



^Jriririr + 1) 



/ 1 \ 



Introducing the notation 



/ mi/ti - 1 \ 



can be written as 



X 



al + {rir + 1) 



+ 



i=2 



a2 + ( a^^ 



(Bl) 



+ 



cr,1 



From Eq. ()B1|) we have y ■ ui = D(Lf)/A/n7H-T, y ■ Uj = {j = 2, ■ ■ ■ ,nr), y ■ Un^+i 
— \Jnr/ {rir + l)-D(i^f), and we finally get 



< + [rir + 1) 
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APPENDIX C: DERIVATION OF THE SYSTEMATIC LIMIT IN THE KASKA 
PLAN 



It is easy to show that the diagonahzed matrix out of the covariance matrix V (P7j) is 

dmg{crl,al + A+,crl + A^), (CI) 

where 



2 ^ 1/2 

r)N4 



The corresponding eigenvectors are 
1 



u 



■(1) 



-4 
-3 



-A+ + 46r] 



-A_ + 4:6r] 



u 



(2) = A/; -A+ + 53r7 , u^^') = U- -A- + SSr^ 
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-A+ + 49r7 



-A_ + 49r] 



where 



0- 



and A/± are the normahzation constants. Hence we get 



X 



ufD{Li 



~ sin^2^ 



74 a2 



49 [DiLi)f 



a2 + A- 



74 a2 

In the actual KASKA case, using the reference values (0) and (0), from numerical calcu- 
lations we obtain the diagonalized covariance matrix 

diag(Ai, A2, A3) = (T^diag(1.044, 4.651, 90.878), 

and the corresponding three eigenvectors 



u 



(1) 



u 



(1) 

2 

(1) 



-0.4975 \ 
-0.3114 , 
0.8097 / 



0.6503 \ 
-0.7516 , M^^) 
0.1105 / 



0.5741 \ 
0.5815 . 
0.5764 J 



can be written as 



Ai 



A. 



A. 



Ai 
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where y in this case is given by 




- sin^ 29 



^ ''1 

a=l 
a=l 
\ a=l 



a3j 



= -sin^2e L)(L2 



D{L, 



J 



using the definition of D{L) in Eq. Q and D{Lj) in Eq. (jHH). Hence we get 



sin^ 26 



X 



IMal 

from which Eq. (j3U|) follows. 



(1) 
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Reactors (a) 


Power/ GWth 


T 1 






1 


3.293 


482 


1663 


1309 


2 


3.293 


401 


1504 


1224 


3 


3.293 


458 


1374 


1233 


4 


3.293 


524 


1149 


1169 


5 


3.293 


1552 


371 


1484 


6 


3.926 


1419 


333 


1397 


7 


3.926 


1280 


340 


1306 



TABLE I: The Power of the the reactors in GWth ^-nd the distance Laj in meters from the three 
detectors (j=l(near), 2(near), 3(far)) to each reactor (a = I,-- - ,7). The powers of reactors are 
hsted in |http : // www. tep co ■ co . j p /kk-np / index-j . ht ml[ The positions of the reactors were read from 
a map and are subject to a few meters of inaccuracy. However, such inaccuracy hardly affects the 
estimation of the systematic errors in the text. 



Reactors (a) 


tal/tl 


ta2/'t2 




1 


0.208 


0.012 


0.133 


2 


0.301 


0.015 


0.152 


3 


0.231 


0.017 


0.149 


4 


0.176 


0.025 


0.166 


5 


0.020 


0.239 


0.103 


6 


0.029 


0.353 


0.139 


7 


0.035 


0.339 


0.159 



TABLE II: The fraction taj/tj of the yields at each detector j from the reactor a in the KASKA 
plan. 
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L (b) 




(d) 



® 



(e) 




(c) 



® 





(f) 




® 







FIG. 1: Configurations of tlie experiments: (a) reactors + one detector, (b) reactors + 
{rir + 1) detectors, (c) The KASKA plan with 7 reactors + 3 detectors, (d) The ideal limit of 
the KASKA plan, (e) One reactor + 2 detectors to be compared with (c). (f) The hypothetical 
KASKA plan with the near detectors in the wrong locations to be compared with (c). 
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0.1 



CM 



0.01 



0) 

E 

< 



0.001 



0.0001 




90% CL 
( 1 d.O.f. ) 



actual KK 20t yr 
1r+ 2d 20t yr — 



0.01 



0.1 



1 



sin22ei3 



FIG. 2: The comparison of the sensitivity to sin^ 2^13 for the actual KASKA plan (the solid line; 
with the configuration depicted in Fig^c)) and a hypothetical case (the dashed line; with the 
configuration depicted in Fig^e)). The statistical errors as well as all the systematic errors are 
taken into account. 
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0.001 0.01 0.1 1 



sin220.,3 

FIG. 3: The comparison of the locations of the near detectors in the KASKA plan. The solid 
(dashed) line is for the actual (hypothetical) KASKA plan with the location depicted in Fig^c) 
(Fig^f)), respectively. 
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0.016 

0.017 

0.018 

0.019 

0.020 

X optimized 
^ planned 




— 1km 



— 1 2km 



w/o Stat, error |Am3^ |=2.5x10"W 



FIG. 4: The contour plots of the systematic hmit on sin^ 2^13 in the KASKA experiment. The 
optimized and currently planned positions of the detectors are also depicted. When the contour for 
each detector is plotted, it is assumed that other detectors are located in the optimized positions. 
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0.020 

0.021 

0.022 

— 0.023 

0.024 

X optimized 
^ planned 




— 1km 




— 1 2km 



20ton»yr |Am3j=2.5x10"^eV^ 



FIG. 5: The same contour plot of the sensitivity to sin^ 2^13 as Fig.lUwith the data size of 20ton-yr. 
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